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Abstract 

The in-plane polarizational stopping power of heavy-ion diclusters in a 2D strongly coupled 
electron liquid is studied. Analytical expressions for the stopping power of both fast and slow 
projectiles are derived. To go beyond the random-phase approximation we make use of the in- 
verse dielectric function obtained by means of the method of moments and some recent analytical 
expressions for the static local-field correction factor. 
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I. INTRODUCTION 



Investigation of static and dynamic properties of (quasi-) two- dimensional systems of 
charged particles has become an important research area in condensed matter as well as 
in physics of strongly coupled many-body systems. 

Nowadays, it is possible to find a number of examples of experimental realizations of sys- 
tems consisting of electrons which are confined within a two-dimensional configuration. The 
most relevant example in the classical statistical regime could be that of electrons trapped 
in the liquid helium surface.- But a major number of those experimental realizations corre- 
spond to quantal electronic systems^ appearing in the fabrication of nanoelectronic devices, 
such as, e.g., the semiconductor-insulator junctions of semiconductor heterostructures. 

On the other hand, the stopping power permits to characterize the interaction of charged 
particles with matter and, therefore, it has become a useful diagnostic tool for experiments 
with polarizable media. Although, there are several linear and non-linear physical phenom- 
ena contributing to the loss of energy of moving charged particles interacting with condensed 
matter,-'^'^'^'^'^ here we focuse on the part of the energy losses due to the medium polar- 
ization. In this physical picture, the loss of energy of the external projectile is caused by the 
creation of a dipole between the moving projectile and the center of the cloud of induced 
charge which intends to screen the former .= Following the work of Lindhard,^^ the calcu- 
lation of the polarizational stopping power is usually related to the imaginary part of the 
medium inverse dielectric function. 

In the sequel, we consider a two-dimensional electron gas (2DEG) consisting of n electrons 
per unit area immersed in a uniform and rigid neutralizing background of positive charges. 
We will assume electrons to interact via a 3D Coulomb potential inversely proportional to 
the distance. 

In the completely degenerate case, such a system can be characterized by a single param- 
eter, which is proportional to the Wigner-Seitz radius, a, i.e., the Brueckner parameter 



where is the Bohr radius. 

In addition, this Brueckner parameter can be understood as a coupling parameter, i.e., as 
a ratio between the characteristic Coulomb interaction energy, e^/a, and the characteristic 



a 

as 
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kinetic energy which is the Fermi energy, 

1 /2 

kp = {2Txn) and m being the Fermi wavenumber and the electron mass. 

In a finite-temperature system, the characteristic kinetic energy is proportional to the 
temperature in energy units, (3~^, yielding the coupling parameter 

In this case, an additional parameter is needed to describe the electron fiuid, the one which 
quantifies its degeneracy, 

D = pEp. (4) 

The strongly coupled regime is usually defined by a value of the coupling parameter higher 
than unity. Under such conditions, the mean-field theories, such as the random-phase approx- 
imation (RPA), are unable to account for interparticle correlational effects and, generally 
speaking, fail to describe the static and dynamic properties of these systems properly. Here 
the Wigner crystallization^*^ constitutes a natural limitation for our approach. 

The in-plane polarizational stopping power of a two-dimensional Coulomb system was 
considered by Bret and Deutsch, who obtained an expression relating this quantity, for 
the case of a general extended projectile, with the loss function of the system.— Perhaps, 
the most outstanding feature of this problem is that the single-ion stopping power in the 
fast-projectile limiting case exhibits a characteristic asymptotic form which is inversely pro- 
portional to the projectile velocity,^* in contrast to the well-known BBL formula applicable 
in 3D fluids." This result has been conflrmed by other authors by means of the linear 
response theory, the scattering theory,-i^ii^ and the harmonic oscillator model.— 

In particular, for a dicluster-like charged particle distribution, the expression of the in- 
plane polarizational stopping power reduces to^^ 



dE ( dE\ ( dE\ 



- 



dx V d^Jnncorr \ dx 



oo kv 

(-^] = — [Z! + Zl) [dk [ du ^ Im (--^] , (6) 
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, , — 2Z1Z2 / dk / duj ^ ' =lm , (7) 

V \ / 

where the first term accounts for the stopping power of two point-hke uncorrelated external 
projectiles with charge numbers Zi and Z2, respectively, whereas ([7]) stands for the correlated 
contribution. We assume here that the two-point-like projectiles are randomly oriented (see 



Refs. 
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In addition, in Ref. [ij, the 2D RPA dielectric function was derived and the asymptotic 
forms for the polarizational stopping power corresponding to the limiting cases of fast and 
slow single point-like projectiles and diclusters were obtained." Later, in the case of single 
point-like projectiles,— an interpolation formula for the local-field correction factor for a zero- 
temperature 2DEG obtained in Ref. 2l| was used to account for the electronic correlations. 

Here, our main aim is two-fold: first, by means of the method of moments we show 
analytically that the fast dicluster stopping power asymptote of a 2DEG remains unaffected 



by correlational effects, like it was done in Ref. l22| for a single projectile. Second, we revisit 



the low-velocity asymptote at zero-temperature, but now including both the single projectile 
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and the dicluster cases and applying a more recent interpolation formula found in Ref. 
which reproduces correctly the short-range effects,— and, finally, we also consider this 
low-velocity asymptote in a correlated high-temperature system.— 

II. HIGH- VELOCITY ASYMPTOTIC FORMS 
A. Dielectric formalism 

To study the fast-projectile limiting case, we make use of the inverse dielectric function 
derived by means of the method of moments and the Nevanlinna formula.— 1^ The same 
approach has been applied in the treatment of the 3D and 2D electron gas high-velocity 
stopping power asymptote,— 1^ in the case of single-ion projectiles. In particular, in the 2D 
case we can write:— 

(k Z) = l+ ^2dW [z + qik,z)] .g. 

^ ' ^ ^ z[z^-iul{k)] + q{k,z)[z^-Luj{k)y ^ ^ 

where 

ujI (k) = C2 (k) /Co (k) = (1 - e-' (k, 0))-' , (9) 
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ul (k) = C4 (k) /C2 (k) = uoIj, {k) {l + K{k) + L {k)) , (10) 

/oo 
uj''C{k,uj)duj, z/ = 0,2,4, (11) 
-00 

being the three non- vanishing and non- diverging frequency power moments (sum-rules) sat- 
isfied by the loss function C{k,uj) = a;^^Im {—e^^ (/c, a;)):— 

Co(A;) = 7r(l-£-i(A;,0)), (12) 

C2 (k) = TT ^^^^ k = TTCJon, (13) 

m 

where uj2d = ^2D (k) is the 2D plasma frequency, while 

C4 (A;) = Toj^^ {1 + K{k) + L (k)) , (14) 

with 



K (k) = l^^k + -^k' (15) 
27me^ Srcne'^m 



being the kinetic contribution to the fourth sum-rule (-Efcm is the average kinetic energy per 
electron) and 

q(7^o,^k) ^ 
being the correlational contribution. 

The parameter function q {k, z) is an analytic (in the upper half-plane) Nevanlinna-class 
function with some specific mathematical properties.— From a phenomenological view- 
point the easiest way to choose this function is to put q{k,z) = iO~^. This is equivalent 
to assume the existence of a perfectly defined single collective excitation, which acts as the 
main energy tranfer channel from the projectile to the plasma. Under this Ansatz we should 
write the loss function as 



C {k, u) = TT 



—^d {uj + UJ2) H 2 ^ (^) + TT^*^ - ^2) 

Z,UJ2 



(17) 



which resembles the Feynman approximation for the dynamic structure factor.— This 
loss function (fT7|) not only satisfies all three sum- rules (11211141) . but also describes a collec- 
tive excitation of frequency uj2{k) which incorporates correlational effects (by means of the 
contribution (fT6|) ) beyond the RPA.— 
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B. Stopping power 



Physically, the applicability of the canonical form (fT71) of the loss function relies on the 
fact that for the fast-projectile asymptote the main energy transfer mechanism from the 
external ion consists in the creation of plasmons in the dielectric medium. This assumption 
permits to handle this stopping power asymptote analytically even beyond the RPA. 

In fact, by applying the loss function f|T71) . it has been established that the fast-projectile 
asymptote of the polarizational stopping power for point-like projectiles 



dx J r,an V ' 

/ "»icorr * is 



(18) 



which coincides with the result obtained by Bret and Deutscb^*^ in the random-phase 
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are of a more general nature and are 



approximation. Indeed, the arguments used in Ref. 
straightforwardly applicable to an external projectile with a general extended distribution 
either in the 3D case or in the 2D one.— 

Following the derivation of Ref. 3, in the random-phase approximation we can write the 
collective excitation frequency appearing in (fTTl) as 



Ulr, = LU. 



2D 



k + - ^fc^ 



(19) 



27™e^ Snne^m 

where we neglect correlational contributions in (I10|17p . Typically, this approximation is 
applicable in the weakly 2DEG. 

If we account only for the leading terms of the long- and short-wavelength asymptotes, 
we get 



;2 ^ uj^j, 



1 + :^ 



(20) 



This expression together with f|T71) can be introduced into ([7]) to give 



_dE\ ^ 2v/2(2ZiZ2)e^ _^Joi2kpRz)_^^ 

where 6 = {2kpaB)^^, b=vp/v. In addition, -2mm, Zmaxi z^in < Zmax, stand for the two 
positive roots of the polynomial equation — b~'^z + 5 = 0,— which has one negative and 

two positive roots, Zmmi Zmaxj Zmin ^ Zjnaxj if 
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b < 5-2/3. (22) 

The asymptotic expansion of the correlated contribution of the dicluster stopping power 
in 2D (EH) gives 



which coincides with the result obtained in Ref. 



19 and where 



- 2kF^ ~ 2mv ^^^^ 

^ Vp 

is the fast-projectile coagulation distance.—"^ In the expansion (1231) we have assumed the 
ratio r/b, with r = kpR and 6 = vp/v to be finite. But it is clear that when the intercluster 
distance becomes larger than the Fermi radius kp^, the asymptotic form (123|) is proportional 
to v~^. 

In case of a strongly coupled electron fluid, one needs to go beyond the RPA and account 
for interparticle correlations, i.e., to consider the expression for the collective excitation fre- 
quency UJ2 {k) steming from ffTOj) . Here we can again make use of an interpolation formula be- 
tween the long- and short- wavelength asymptotes, {k | 0) and {k f 00), respectively:— 



1 + — ^k + ^k + —k^ 



^2 (k) ^ Uj'^D 

Ec being the correlation energy per particle. 

If we introduce expression f|T7|) . but now with fl25l) . into ([7]), then 

dE 2V2{2ZiZ2)e^ T™"^ b^Jo{2kFRz) 



(25) 



b^Jo{2kFRz) 
A™.,:„ Jl - eb^ - 6bH-^ - b^z^ ' 



dx Tsal - _ 6bH-^ - b^z^ 

with 

which accounts for the contribution of the average kinetic and correlation energies per elec- 
tron. 

Obviously, this integral can be recast into 
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with 



and where condition (1221) is replaced by 



b < , ^ (30) 

Clearly, the leading term of the asymptotic expansion of ( l28l) coincides with ( l23l) . Therefore, 
correlational effects beyond this approach, which are accounted for in the canonical form 
(fT7|) by means of (fT6!) . do not affect the asymptotic behavior of the polarizational stopping 
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power obtained in Refs. 

The expressions ( IT8|23l) have no finite limit when h — > 0.— Notice that the contribution 
to the dispersion law steming from the electron-hole excitation in ( fT5l) . which is propor- 
tional to fc^, is the main factor responsible for the convergence of the integrals ( |6|71) in the 
fast projectile limiting case.— Then, even for a high-temperature electron liquid, the 
minimum Coulomb impact parameter which is usually introduced as an upper cut-off in the 
integration over the wavenumber, 

/i being the reduced mass, exceeds the de Broglie wavelength and, hence, must be replaced 
by the latter,— i.e., 

However, the inclusion of the electron-hole excitation term guarantees the convergence of 
the above-mentioned integrals without introducing the cut-off. 



III. LOW- VELOCITY ASYMPTOTIC FORMS 



Apparently, the stopping of low-velocity projectiles due to linear polarizational effects 
was first studied by Fermi and Teller,— who found a characteristic linear dependence of the 
stopping power on the projectile velocity. This result seems to be valid for any degree of 
degeneracy of the electron fluid and any coupling regime. In fact, this is just a consequence 
of the heavy projectile approximation, M — ^ oo}^ It must be pointed out that, since we are 
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assuming the trajectory of the projectile to be a straight hne, its kinetic energy must always 
be much greater than the kinetic energy of the electrons, or, equivalenty, 

/ m 
V M 

Thus, when we refer to the slow projectile limit we must assume the condition 

/ m V 
V M vf 

to hold.i^ 



A. Dielectric formalism 

In the slow-projectile limiting case, contrary to the previous one, we expect the velocity 
of the projectile to be too small to create a plasmon in the electron liquid. Thus, now the 
energy transfer to the plasma is essentially determined by the low-frequency range of the 
fluctuation spectrum. In addition, in spite of the efforts to find an adequate expression 
for the parameter function q{k,uj) appearing in ([8]), capable to describe this range of the 
spectrum^ adequately, here we make use of the general static local-field-corrected dielectric 
function, 

e{k,uj) = l + (P{k)U{k,uj) = l + - .... Tjr^, (33 

which permits to apply recent numerical results on the static properties of the electron 
gas obtained by numerical simulations. In the previous expression (j){k) = 2'Ke^/k is the 
two-dimensional Coulomb potential, 11 (A;, uj) is the polarization function of the system, 
n_RPA (fc, oj) its form within the RPA, and G (k) is the static local-field correction function. 
Further, from expression fl33p we can derive the loss function: 

aku;) = ^ <^if')^nPAik,u;) 

^'^ uj[l + 4>{k)H{k)U'^p^{k,u;)f + [^{k)H{k)U'i,p^{k,u;)f 
where we define the function H {k) = 1—G {k) and Hrpa {k, u) = H'ppA {k, Lj)+ill'^p^ {k, u). 

Obviously, to recover the RPA we must put H {k) = 1. 

B. Stopping power 



Next, we introduce this last expression into the general equations ([6]) and ([7]). In or- 
der to simplify the corresponding integrals, we can make use of two common approximations: 

9 



(i) since we are dealing with slow projectiles, v/vp ^ 1, we substitute the denominator of the 
loss function by its static limiting form {u = 0) and recall that H'^p^ [k, 0) = 0; (ii) we can 
approximate in the numerator the imaginary part of the polarization function H'^p^ {k, u) 



1. Zero-temperature system 

In the zero-temperature limiting case, the slow-projectile expression corresponding to the 
point-like contribution ([H]), can be recast as 



_dE\ 



2V2{Zf + Zi)e^ V 



s^B 



dz 







which corresponds to that derived in Ref. 
contribution (171) we have 
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(35) 



for a point-like projectile, whereas for the 



dE\ 
dx J 



2V2{2ZiZ2)e^ 



Vf 



dz z'^Jq {2kFRz) 



z + ^,H{z) 



2 ■ 



(36) 



These two expressions generalize those derived by Bret and Deutsch for slow projectiles.— 
Again, these are recovered under the RPA, i.e., with H (k) = 1. 

On the other hand, non-linear (non-perturbative) stopping power effects beyond our (lin- 
ear) polarizational approach might be important in the slow-projectile limiting case. These 
effects have been treated using the non-perturbative scattering theory in a partial-wave 
representation, both for single heavy ions->^"^i^ and diclusters.- In particular, for randomly 
oriented slow diclusters of similar 'atoms', in Ref. 19| it was shown that the correlated contri- 
bution to the differential cross section, appears via a Bessel function factor similar to that 

of (IHSD. 

To go beyond the RPA we need an analytical expression for the local-field correction 
factor, Gik), of a paramagnetic two-dimensional degenerate electron gas. For instance. 



f. y 



m Ref. 1161 an interpolation formula obtained in Ref. 



21 



was used. In additon, a more 



recent interpolation formula can be found in Ref. 



23l . Indeed, this latter expression not only 



accounts for the well-known long- wavelength asymptotic behavior of the local-field correction 
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factor, but also reproduces correctly the short-range effects.— In particular, 



G(A;iO)^AA, (37) 



and the constant A is determined via the compressibility sum-rule to give^ 

^=^fl--)^ (38) 



V2rs V K 

where kq = 7rrf/2 is the compressibility of the ideal electron gas (in the units of a%/Ry), 
whereas k is the compressibility of the interacting gas. For T = 0, we have 

^ = 1 _ + !l ( d'Ecjrs) _ 1 dE,ir,) \ _ 

K TT 8 \ drg"^ Tg drg J 

Here, for the correlation energy per particle, Ec, one might use, for instance, the expres- 
sion obtained by Rapisarda and Senatore^ by fitting of data of the diffusion-Monte-Carlo 
simulations or the one fitted from the quantum-Monte-Carlo simulation results in Ref. 3^ : 



e2 



-0.1925 + (0.0863136r, + 0.0572384r2 + 0.00363rf) (40) 



X In 



1 + (l.0022r, - 0.02069rf 2 + 0.33997r2 + 0.01747r,^)" 



which has been obtained for 1 < < 40 (see also Ref. l35| and references therein for a 
further discussion). 

On the other hand, for the short-wavelength asymptotic form we have^^i^^ 

G (A; t oo) ~ C-^ + 5, (41) 

where C depends linearly on the difference in kinetic energy between the interacting and 
the ideal electron gas,— 

^^-57i^<^'^=<^'>'' '''' 

In addition, B = 1 — g (0), where g (0) is the value of the pair-correlation function at the 
origin. For the latter we can make use of the interpolation formulai^^ 



V2_ 

l + 1.372r, + 0.0830r2 



^(0)= . , .o.o.^:..oo... - (43) 



23l for the static local-field correction 
. I33 valid in the 



Although the interpolation formula obtained in Ref. 
factor was explictly derived by using the correlation energy formula of Ref. 
range 0.1 < < 10, here we also use this expression with the interpolation formula (HOj) for 
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the correlation energy obtained in Ref. |3J. The results for the static local-field correction 
are compared in Fig. [H Notice that the expressions obtained by means of the interpolation 
formula of Ref. 



23 



37 



show the shift of the peaks beyond 2kp, as decribed in Ref. 
In Fig. [2]we plot the ratio between the low-velocity asymptotic form of the stopping power 
of a single projectile, i.e., the sum of (153]) and with _R = 0, and its RPA counterpart. 



c 



_dE\ ( _dE\ 
dx J \ dx J 



-1 



> 1, 



RPA 



and study the difference obtained using the interpolation formula of Ref. 
expression for the correlation energy of Ref. 



23 



33 



or of Ref. 



(44) 
(with the 
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as it was done in Ref. 



16 



341 . respectively) or that of Ref. 



As one can see, the electronic correlations enhance the stopping power, this increase 
becoming more pronounced for lower electron densities (high coupling region), as it was 
outlined previously.— 

Moreover, although we only need to integrate in fl35ll36p until the value oik = 2kp,we see 
that there is a certain quantitative discrepancy due to the usage of the expressions for the 



correlation energy Ec (r^) from Ref. 



33 



or from Ref. 



34. In both cases, it is also noticeable 



that the denominator in the integrand of fl35|36l) vanishes for the values of > 10, which is 
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indeed beyond the fitting range of G (k) in Refs. |21 

It is clear from (136|) that in the slow-projectile limiting case the natural scale for the dis- 
tance between the ions in the dicluster is the coagulation distance in a completely degenerate 
system,—"^ 



R' 



2 



as 



(45) 



2kF 2V2 

In Fig. [3] we compare the ratio of the correlated low- velocity asymptotic form (!36|) to the 
uncorrelated one fl35|) of a dicluster distribution as a function of the distance between the 
ions in atomic units, for different values of the Brueckner parameter, in both the RPA and 
the beyond-the-RPA cases, and assuming that both ions have the same charge number. 

In particular, an increase of the value of the Brueckner parameter makes the Friedel- 
like oscillatory pattern of ( l36i) to exhibit a lower damping rate (Fig. [3] (a)) due to the 
proportionalty between the Brueckner parameter and the coagulation distance R'^ f H5l) . 

In addition, the comparison of the above-mentioned ratio of the correlated low-velocity 
asymptotic form fl36l) to the uncorrelated one fl35l) in the RPA and beyond the RPA, might 
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provide a quantitative insight into the dependence of the screening distance on the value of 
the Brueckner parameter (Fig. [3] (b)). 

The 2D Thomas-Fermi wavelength, Xtf = 0^/2,— represents a typical screening distance 
in the random-phase approximation. In the region where the RPA is valid, ^ 1, the 
Wigner-Seitz radius is smaller than as- This means that the total number of particles which 
are contained inside the disk of radius Xtf, 

N = T^Xl^n = -L, (46) 

s 

is large enough to screen the potential created by the external ion effectively. On the other 
hand, beyond the domain of applicability of the RPA, > 1, the number of particles 
inside the disk of radius Xtf is smaller than unity and, then, the Thomas-Fermi wavelength 
becomes an underestimate of the typical screening distance. Thus, we observe that the 
oscillations exhibited at = 10 in the RPA case in Fig. [3] (b) become even less damped 
when the local-field correction factor G{k), which permits to describe the behavior of the 
system at higher values of the Brueckner parameter more accurately, is accounted for. 



2. High-temperature system 



In the high-temperature limiting case, expression ([6]) with reduces to 



dx J 



oo 

a/tt (Z^ + Z2) pvki) J dx- 



x'^ exp (— x^) 



X 



V2 ^ 



1 2 ' 



(47) 



where 

as 



2'ime^f3 is the 2D Debye wavenumber, p = me'^/h, whereas ([7]) can be written 



_dE\ 
dx) ^ 



2y/TTZiZ2pvk£, I dx 



exp 



') Jo [Rx) 



X -\- 



V2 



H 



2 ' 



(48) 



with R = 2R/X, X = (h^P /2mY^'^ being the thermal de Broglie wavelength. Again, if we 
put H (k) = 1 we recover the result obtained by Bret and Deutsch for the limiting form of 
the stopping power of slow projectiles^*^ of a high-temperature system. 

Once more, Eqs. fl47|48l) have no finite limit as h —>■ 0, although now the previous cut- 
off argument is not applicable. Nevertheless, one migth argue that to avoid the Coulomb 
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collapse on a slow projectile, one has to take into account the effects of electronic diffraction 
essentially related to the finiteness of h. Notice that this peculiarity, in both the fast- 
projectile limiting case and the slow-projectile one, is shared by the 3D problem as well,— 
and, therefore, seems to be specific for the calculation of the polarizational contribution to 
the stopping power. 

Since the expressions for the local-field correction factor employed before are only appli- 
cable to zero-temperature systems, to go beyond the RPA we consider here an interpolation 
formula based on the known asymptotes that G (k) must fulfill.— 

First, the long- wavelength behavior is described by expression ( 1371) . being the value of the 
constant A again dictated by the compressibility sum-rule. For a high-temperature system, 
D < 1, we have^ 



where for the correlation energy per electron of a classical electron gas we can use the MC 
formula of Totsuji^ 

pEc (F) = -1.12F + O.riF^/^ - 0.38, (50) 

which is valid for < F < 50. 

For the long-wavelength behavior of the local-field correction factor in the high- 
temperature system, one expects the Holas contribution^^ to vanish. Indeed, in the classical 
limiting case the kinetic energy of an interacting electron system coincides with that of a 
non-interacting one. Hence, we can write 

G(A;r 00)^1-^7(0), (51) 

where we can even assume that the electron-electron pair-correlation function vanishes at the 
origin in the high-temperature limiting case. With this in mind, the interpolated expression 
for the local-field correction factor can be cast as 

G{z) = ^^, (52) 

where z = k/2kp. This high-temperature static local-field correction is used to plot Figs. H] 
and [5] for a dicluster projectile with Zi = Z2. 

In Fig. m we display the ratio (jBl) for a high-temperature strongly coupled 2D system at 
different values of the degeneracy parameter D. As in the completely degenerate case before, 
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one can see that ( (F; D) > 1, and the departure from the RPA becomes more notorious as 
the couphng parameter increases, although the former is shghtly smaller at lower values of 
the degeneracy parameter than the one found in the zero-temperature system, i.e., there is 
a consistent tendency of enlargement of the ratio C (F; D) as the degeneracy increases. 

In the classical limiting case, -D ^ 1, the natural length to measure the distance between 
the cluster ions (see Eq. (136|) ) is substituted by the thermal de Broglie wavelength which is 
proportional to the coagulation distance in the slow-projectile limiting ca.se,— ^ 

K = - = (53) 

since it represents the minimum distance which might be resolved in the electron liquid due 
to the uncertainty principle. Here the typical screening distance in a weakly coupled plasma, 
e.g., the Debye radius, is also much larger than the thermal de Broglie wavelength. 

Further, due to the relationship expressed in (1531) . the behavior of the correlational term 
f HHj) is similar to that of the zero-temperature system, i.e., its damping rate decreases as the 
value of the coupling parameter increases, although this tendency appears to be obviously 
compensated at higher values of the degeneracy (Fig. [S]). It is also noticeable that when the 
degeneracy rises the contribution psj) becomes negative and there is a tendency to reproduce 
the Friedel-like oscillations which characterize the completely degenerate system. 

As in the previous case, the evaluation of the dependency of the screening distance on 
the value of the coupling in the RPA and with a local-field-corrected dielectric function, 
clearly indicates that the quantitative difference found previously remains applicable at 
high temperatures. Indeed, under these conditions the number of particles inside the disk 
of radius A^), 

^ = ^^D^ = 4fl' (54) 
is insufficient to screen the disturbance induced by the ion potential, and the typical Debye 
radius becomes again an underestimate of the screening distance as the coupling increases. 



IV. CONCLUSIONS 



Correlational contributions to the in-plane polarizational stopping power of heavy-ion 
diclusters by 2D strongly coupled electron fluids have been assessed. The limiting forms for 
fast- and slow-projectiles and for both low- and high-temperature cases have been considered. 
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In the case of a high-velocity projectile we have used a dielectric formalism based on the 
employment of the canonical solution of the truncated problem of moments for the loss 
function, and have established that the fast-dicluster stopping power asymptote is unaffected 
by the interparticle correlations. On the other hand, for slow projectiles we have compared 
several fitted expressions for the static local-field factor. In particular, an ad-hoc interpolated 
formula has been constructed for the classical system. In this slow-projectile asymptote the 
correlational effects are shown to be quantitatively important, as it was previously outlined. 
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Figure Captions 



Fig. 1. The fitted expression for tfie static local-field correction factor of: (a) Ref. 

n . n 



(b) Ref. 35 with E^. (r^) from Ref. 
r„ = 0.1; dash-dot: r„ = 1: solid: r 



33 



21 



; (c) Ref. l35| with Ec (r^) from Ref. '341 Dash-dot-dot: 
5; dashed: r, = 10. 



Fig. 2. Ratio between the low-velocity asymptotic expression for the stopping power 
and its RPA counterpart (14^ in the zero-temperature limit and for a single projectile, i.e., 
R = 0. The fitted expression for the static local-field correction factor of Ref. |35| with 



Ec (r J from Ref. l33| (dashed) and from Ref. [3^ (solid) are compared to the expression of 



Ref. 



21 



(dash-dot). 



Fig 3. Ratio between the correlated low-velocity asymptotic expression ( !36|) for the 
dicluster stopping power and its uncorrelated counterpart ( !35l) in the zero-temperature 
limiting case with Zi = Z2. (a) RPA for = 2 (solid) and = 10 (dash-dot); (b) 
Comparison of different fitted expressions for the static local-field correction factor for 



10: RPA, i.e., = 0, (solid), G{k) o 



dot-dot), G{k) of Ref. 



35 



with ErirS] from Ref. 



Ref. |35| with Ec (r^) from Ref. 33| (dash- 
34 (dash-dot), and G{k) of Ref. I2I (dashed). 



Fig. 4. Ratio between the correlated low-velocity asymptotic expression for the stopping 
power and its RPA counterpart (jBl) in the high-temperature limiting cases and for a single 
projectile, i.e., = 0. The interpolated expression for the static local- field correction 
factor (!52|) is used. The values of the degeneracy parameter, are: 10"^ (solid), 2 x 10^^^ 
(dash-dot), and 10~^ (dashed). 



Fig. 5. Ratio between the correlated low-velocity asymptotic expression fHHl) for the 
dicluster stopping power and its uncorrelated counterpart (H7j) in the high-temperature 
limiting case with Zi = Z2 for D = 10^^ (a) and D = 10^^ (b). The values of the coupling 
parameter are: F = 2 (with G{z) = 0, (solid), or G{z) from fl52|) . (dash-dot)) and F = 10 
(with G{z) = 0, (dash-dot-dot), or G{z) from ([52D, (dashed)). 
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